Quantum rattling of molecular hydrogen in clathrate hydrate nanocavities 
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We have performed high-resolution inelastic neutron scattering studies on three samples of hydro- 
genated tetrahydrofuran-water clathrates, containing either H2 at different para/ortho concentrtion, 
or HD. By a refined analysis of the data, we are able to assign the spectral bands to rotational and 
center-of-mass translational transitions of either para- or ortho-H2. The H2 molecule rotates almost 
freely, while performing a translational motion (rattling) in the nanometric-size cage, resulting a 
paradigmatic example of quantum dynamics in a non-harmonic potential well. Both the H2 rota- 
tional transition and the fundamental of the rattling transition split into triplets, having different 
separation. The splitting is a consequence of a substantial anisotropy of the environment with 
respect to the orientation of the molecule in the cage, in the first case, or with respect to the center- 
of-mass position inside the cage, in the second case. The values of the transition frequencies and 
band intensities have been quantitatively related to the details of the interaction potential between 
H2 and the water molecules, with a very good agreement. 

PACS numbers: 78.70.Nx, 82.75.-z, 63.20.Pw 



In the presence of a low amount of other molecular sub- 
stances, frozen water can form several crystalline com- 
pounds, different from common ice, known as clathrate 
hydrates The hydrogen-bonded water lattice can 

host different guest molecules in cages of various sizes 
and geometries. The deeply non-harmonic motion of the 
guest molecules differs from all other normal modes of 
the crystal, and is responsible for the anomalous thermal 
conductivity of these materials Q • Besides being present 
in nature, these compounds can be prepared in labora- 
tory and have been recently proposed as effective, safe, 
and economical materials for hydrogen storage 0, 0, ■ 
Binary hydrogen clathrate hydrates (i.e. made of H2O 
and H2, only) require ~ 2000 bar of pressure to be pro- 
duced at T ~ 273 K. However, it has recently been shown 
that the ternary compound with tetrahydrofuran (THF) , 
can store significant amounts of hydrogen with a much 
lower formation pressure 0]. An understanding of the 
interaction of the H2 molecule with the host material is 
fundamental for the rational design of clathrates as hy- 
drogen storage materials. This can be accessed by an 
inelastic neutron scattering (INS) study of the dynamics 
of the molecule trapped in the cage, as presented in this 
work. 

The crystal structure of the II2O-TIIF-H2 clathrate is 
cubic (sll structure) with 136 II2O molecules in the unit 
cell, giving rise to sixteen (small) dodecahedral cages and 
eight (large) hexakaidecahedral cages 0, The THF 
molecules occupy the larger cages, and no more than one 
hydrogen molecule is hosted in the small cages [7| . There 
is an intrinsic advantage in using INS with hydrogenous 
materials due to the incoherent scattering cross-section 
of the proton, which is almost two orders of magnitude 
greater than in any other nucleus, thus allowing for a rel- 
atively simple access to the self dynamics of molecular hy- 



drogen. Hydrogen clathrates for the present experiment 
were produced at ISC~CNR using D2O and completely 
deuterated tetrahydrofuran (TDF) in stoichiometric pro- 
portion (17:1 niol). Hydrogen gas at ~ 800 bar was 
added to the sample by either letting the liquid mixture 
to freeze in the presence of H2, (at T ~ -|-2°C) or al- 
lowing gaseous H2 to diffuse into the ground D2O-TDF 
clathrate (at T ~ — 10°C). Gas-release thermodynamic 
measurements gave results consistent with the hypothesis 
of single H2 occupancy of the totality of the small cages. 
The neutron measurements were carried out on three gas- 
charged samples and one reference sample of pure D2O- 
TDF clathrate. Two samples contained H2 at different 
ortho-para (0-H2 - P-H2)) concentrations, and one sample 
contained HD. Due to the procedure of sample prepara- 
tion, the final 0-H2 concentration of the clathrates ob- 
tained starting with normal-}i2 and almost pure P-H2 
resulted 53 % and 48 %, respectively (o-rich sample and 
p-rich sample). Raman measurements performed at ISC- 
CNR before and after the neutron experiment provided 
an independent determination of the ortho-para ratio in 
our hydrogenated samples. 

For the present experiment, performed at T=20 K on 
the high resolution {Alu/uj < 2.5%) TOSCA spectrome- 
ter at ISIS (U.K.), the interesting energy range is between 
3.5 and 120 meV. Many features due to acoustic, optic, 
and molecular excitations of the lattice were recognizable 
in the spectrum of the reference sample. The spectra 
from the hydrogen-charged clathrates showed, in addi- 
tion, the presence of several narrow and intense bands 
due to the molecular hydrogen dynamics. The hydrogen 
contribution was obtained by subtraction of the reference 
spectrum, which has a low intensity ( 5-;- 15 % of the main 
bands peaks). It was possible to interpret the observed 
lines unambiguously as being due to the molecular ro- 
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tation (E ~ 14.7 meV), to the fundamental transition 
of the quantum ratthng (E ~ 10.0 meV), and to their 
combinations, as discussed below. 

Molecular hydrogen trapped in the clathrate cages ap- 
pears, in general, as a non-equilibrium mixture of 0-H2 
and P-H2 . Thus, we need to consider the neutron scatter- 
ing cross section for the two species, and its dependence 
on the rotational transitions [8| . Neglecting the coherent 
part of the scattering, on account of the overwhelming in- 
coherent scattering length for the proton, the double dif- 
ferential neutron scattering cross section is proportional 
to the self-part of the dynamical structure factor for the 
motion of the center of mass (c.o.m.) SgcitiQ,'^) and is 
written as Q 
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where fci and k{ denote the initial and final neutron 
wavevectors, lu = E/h, and the symbol (g) denotes a con- 
volution product. This expression holds if the hypothesis 
of decoupling of the internal (rotational) motion of the 
molecule from the c.o.m. motion is satisfied. The Dirac 
(^-functions are centered at the energies of the rotational 
transition, hcujj' of the H2 molecule and the intensity 
factor iy(J, J', Q) is a function of the momentum transfer 
Q and depends on the rotational transition J^J' of the 
molecule. The expected spectrum therefore consists of 
the c.o.m. excitations replicated and shifted by the en- 
ergy of any possible rotational excitation of the molecule. 
Since INS is not subject to selection rules, all transitions 
to molecular rotation and c.o.m. vibration states are 
allowed. The intensity factors i'{J, J', Q) have been cal- 
culated assuming a rigid rotor molecule [l3]- At the low 
temperature of the experiment, only the lowest rotational 
states (i.e. J = 0, 1) are populated, and few transitions 
contribute to the spectrum in the observed frequency re- 
gion, namely the rotationally-elastic J = 1 ^ 1 and the 
inelastic J = 1 ^ 2 transition of 0-H2 and the inelastic 
J=0— >1 of P-H2 . We would note that the J=0— *-0 transi- 
tion of the P-H2 molecule, being weighted by the coherent 
scattering length only, does not contribute appreciably to 
the spectrum and, consequently, the observed fundamen- 
tal rattling transition (E ~ 10.0 meV) is mainly due to 
the ortho molecules. 

For unambiguous assignment of all of the spectral fea- 
tures, we compare the spectra of the o-rich and j»-rich 
samples. The intensity ratio of the two bands at ~10 
meV and ~14 meV is lower in the p-rich sample, con- 
firming that the 10 meV band is due to 0-H2 (namely, to 
the combination of the elastic J=l— i'l rotational transi- 
tion with the fundamental rattling transition of the c.o.m. 
motion), and the 14 meV band to P-H2 (namely, to the 
inelastic J=0— ^l rotational transition). By combining 
linearly the p-rich and o-rich spectra, we are able to ob- 
tain the spectra of pure p-il2 and pure 0-H2 . The result is 
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FIG. 1: (Color online) Upper panel: The clathrate spectrum 
(blue line, vertically displaced) can be experimentally decom- 
posed into the sum of the 0-H2 and P-H2 contributions (black 
and red line respectively). The spectrum of P-H2 shows the 
strong J = — > 1 rotational band (a triplet at ~ 14 meV) , 
while the 0-H2 spectrum displays the (split) band around 10 
meV, due to the quantum rattling motion of the molecule in 
the cage. In the lower panel the P-II2 (red) and 0-II2 spectra 
(green) have been redrawn shifted by —AEqi — —14.7 meV 
and by — AiJi2 — —29.0 meV, respectively, to show the co- 
incidence, in frequency and shape, of the rotational-rattling 
combination bands with the rattling fundamental and first 
overtone. 



particularly elucidating and is shown in Fig. 1. The spec- 
trum of 0-H2 (black line) contains contributions from the 
fundamental transition of rattling vibration, and from its 
overtones, while the P-H2 spectrum (red line) contains 
the pure rotational (J=0^1) transition at about 14 meV 
plus combinations of this one with the c.o.m. transitions. 
This can be demonstrated just by shifting the P-H2 spec- 
trum by the energy of the rotational transition (red line 
in the lower panel of Fig. 1) obtaining an almost perfect 
match of the translational bands. The 0-H2 spectrum 
contains also the band due to the J = 1 ^ 2 transition 
(about 29 meV) and the combination of this one with 
the c.o.m. spectrum (shifted spectrum, green line in the 
lower panel of Fig. 1). 

For comparison with the experimental determination, 
we have calculated the spectral intensity of the rattling 
excitations by solving numerically the Schrodinger equa- 
tion for the c.o.m. motion. The potential energy for 
one H2 molecule, as a function of the c.o.m. displace- 
ment R from the center of the dodecahedral cage, has 
been calculated summing the pair H2-D2O potential (as- 
sumed equal to the H2-H2O potential ll|, II4I) over 514 
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FIG. 2: (Color online) Left panel: Isotropic part of the poten- 
tial energy for one H2 molecule as a function of the distance R 
(atomic units) from the center of the dodecahedral cage, and 
calculated energy levels. Right panel: Calculated incoherent 
INS spectra for one H2 molecule (Eq. 1). The spectrum (blue 
line) is the sum of three contributions, where the transition 
lines of the c.o.m. motion (labeled by the final values of A'^, L) 
are shifted by the amount Ajj/ of each rotational transition: 
All = (black line), Aoi = 14.7 meV (red line), A12 = 29.0 
meV (green line). 
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FIG. 3: (Color online) Comparison of the 0-H2 and P-II2 ex- 
perimental spectra with the results of our calculation. We 
have considered the J=l^l and J=1^2 molecular rotational 
transitions for the 0-H2 spectrum and the J=0— >1 transition 
for the P-H2 spectrum. The various bands are assigned. 



molecules of D2O around the center of the cage on a rigid 
lattice, and averaging over H2 molecular orientations. 



The disorder of the water deuterons has been taken into 
account performing several averages over random config- 
urations, all of them respecting the Bernal and Fowler 
ice rules [13] ■ A further average over the direction of R 
gives the isotropic part of the potential. The potential 
energy of the hydrogen molecule as a function of i? = |.R|, 
obtained usin g th e two potential models available from 
literature coincide for the purpose of this analy- 

sis, and is represented in Fig. 2. It is strongly anharmonic 
and rather flat in the center. As a consequence, the cal- 
culated energy levels, (labeled by the principal quantum 
number N — 1,2,3,...), are split (the energy depends 
on the orbital quantum number L — 0,1, 2, ...) and not 
equally spaced. The residual (2L -I- l)-fold degeneracy 
of each level is not removed in this model due to the as- 
sumed isotropy of the potential. The knowledge of the 
wave function for the translation degrees of freedom, and 
the use of equation (1), where all coefficients v{J,J',Q) 
are known, permits the calculation of the intensity of 
transitions, determining the whole neutron energy loss 
spectrum. For each rotational transition JrightarrotuJ' 
of the molecule, we obtain a spectrum given by the rat- 
tling transition energies plus the rotational energy transi- 
tion AEjji. The calculated spectra are shown in Fig. 2, 
where all transitions are represented by Gaussian lines 
having the same (arbitrary) width. In Fig. 3, we com- 
pare separately the 0-H2 and P-H2 experimental spectra 
with the results of our isotropic model, finding a quanti- 
tatively good level of agreement, taking into account that 
this simple model cannot reproduce the splitting of the 
bands. 

The assignment of the bands is confirmed by the re- 
sults obtained for the HD-clathrate, compared to the 
spectrum of the Il2-clathrate in Fig. 4. We also observe 
that the experimental band due to the fundamental of 
the rattling mode (~ 10 meV in the II2 spectrum) is 
split in more than two components, evidencing an ap- 
preciable anisotropy of the c.o.m. potential in the cage. 
In contrast, the splitting of the J=0^1 rotational tran- 
sition into a triplet (~ 14 meV in the H2 spectrum) is 
a consequence of the anisotropy with respect to the ori- 
entation of the II2 molecule. The cage shape, as it re- 
sults from the structural measurements 0, Il4j , is quite 
anisotropic, with the oxygen atoms located at three dif- 
ferent distances from the center. Fitting each band with 
three Voigt lineshapcs (with a Gaussian width fixed by 
the instrumental resolution), we obtain the experimental 
values reported in table 1. Comparing the (average) ex- 
perimental frequencies with the results of our model we 
obtain a very good agreement. In particular, the ratio of 
the rattling frequencies measured for H2 and HD is about 
1.40, and is well reproduced by our calculation. Note that 
this value is intermediate between the harmonic value 
\/ "^HD / ^7iH2 = ^73/2 ~ 1.225 and the value expected 
for a square well potential tohd/"^H2 — 3/2 = 1.5. 

Recently the energy levels of one II2 molecule confined 
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Rattling fundamental 
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8.80 


6.39 


6.59 


1 — 1 
1— 1 


9.94 


8.38 


7.28 


III 


12.53 


9.91 


8.75 


ave. 


9.86 


9.85 8.22 


7.04 7.44 






Rotation 




I 


13.64 


11.32 


10.86 


II 


14.44 


14.66 


11.58 


III 


15.14 


18.83 


12.21 


ave. 


14.41 


14.93 


11.54 



TABLE I: Experimental and calculated transition energies 
(meV). Roman numbers indicate the different components of 
each band, and "ave." means weighted average. 

inside an isolated dodecahedral cage of 20 H2O molecules, 
(similar for dimensions and shape to the one of the THF- 
H2O clathrate), has been calculated, solving numerically 
the five-dimensional Schrodinger equation, and taking 
into account the anisotropy of the potential energy and 
the coupling of rotational and translational motion [ist . 
The authors predicted a splitting into a triplet of both 
the rattling fundamental and the rotational transition, as 
we have experimentally observed (see Table 1). While 
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FIG. 4: (Color online) Rotational and fundamental rattling 
transition for a H2 (lower panel) and HD (upper panel) 
molecule in a THF-clathrate. Each band is fitted with a 
triplet (Voigt shape) whose frequencies are reported in Ta- 
ble 1. 



the calculated splitting for the rattling transition (3.52 
meV) reproduces correctly the experimental one (3.73 
meV), the calculation strongly overestimates the split- 
ting of the rotational transition (7.51 vs 1.50 meV). In 
addition, for the fundamental transition of the rattling 
mode, our calculated result is somewhat more similar to 
the experimental value, probably because we used a much 
larger number of molecules, though using the simplifica- 
tion of an isotropic potential. Then, therefore, the pair 
potential model used in Ref. seems to overestimate 
the actual anisotropic forces on the hydrogen molecule. 

In summary, we have measured the quantum dynam- 
ics of a single H2 molecule in the confined geometry of a 
water clathrate nanocavity. An isotropic model describes 
quantitatively the main features of the motion, were cou- 
pling between rotational and c.o.m. dynamics are weak, 
reproducing frequencies and intensities of the observed 
spectrum. The splitting of the rotational and transla- 
tional bands is a consequence of the anisotropy of the 
environment, which must be modeled with greater accu- 
racy than done until now, to extract direct information 
on the basic interaction between H2 and H2O molecules. 
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